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^ ■ 1. Introduction 



The propagation properties of light beams in nonlocal nonlinear media have attracted sig- 
nificant attention in recent years. There are various interesting properties induced by the 
nonlocality, such as the suppression of collapse [1] and the support of vortex solitons [2] 
and multi-pole solitons [3]. For the special case of the strongly nonlocal nonlinear (SNN) 
media in which the characteristic length of the material response function is much larger 
than the beam width, the propagation equation can be linearized to the well-known Snyder- 
Mitchell model (SMM) [1]. In fact, since Snyder and Mitchell introduced the SMM to inves- 
tigate optical beam propagation in SNN media, various soliton solutions [2|[3|[5]4T0] , such as 



Hermite-Gaussian (HG) [3], Laguerre-Gaussian (LG) [3] and Ince-Gaussian (IG) [6l|9] soli- 
tons, have been predicted theoretically. Some soliton structures and their interaction have 
been observed experimentally in SNN materials such as nematic liquid crystal [Tl][T2]. and 
lead glass [T3] . 

In addition, to the best of our knowledge, the propagation in free space has been inves- 
tigated more thoroughly than the other propagation problems in the field of optics. The 
paraxial diffraction equation has been thorughly investigated. Various types of beam so- 
lutions with different transverse profiles have been obtained in Cartesian, circular cylin- 
drical, and elliptical coordinates (e.g., [TH - ITG] and references therein). These solutions can 
be roughly classified into two types: (i) shape-invariant beams, such as LG, HG, and IG 
beams; and (ii) shape- variant beams, such as higher order elegant-Hermite-Gaussian (EHG), 
elegant-Laguerre-Gaussian (ELG), and elegant-Ince-Gaussian (EIG) beams. The propaga- 
tion of these beams has been investigated in detail and many parameters, such as width, 
divergence, radius of curvature, and quality factor, have been introduced to describe their 
propagation. In summary, the theory of free propagation has been well developed over the 
past decades. 

Compared with that of the free propagation, the problem of SNN propagation is math- 
ematically much more complicated [17], even for the problem of soliton which may be the 
simplest example of the propagation problems in SNN media. However, the structures of the 
HG, LG and IG solitons introduced in the literature are also the modes in free space |18j . 
leading naturally to the question of whether any direct relation exists between free propa- 
gation and propagation in SNN media . 

In fact, if there exists a one-to-one correspondence between free propagation and propaga- 
tion in SNN media, the results of free propagation theory could be applied to the following 
aspects of the study of SNN propagation: (i) The structures of and the existence conditions 
for breathers and solitons could be conveniently described with simple and intuitional physi- 
cal pictures of free propagation in a unified manner, (ii) Whereas the previous investigations 
are mainly focused on solitons and breathers that are shape-invariant upon propagation, the 
shape- variant propagation in SNN media remains unexplored. With an established correspon- 
dence between free and SNN propagation, it would be easy to deal with the propagation of 
an arbitrary field in SNN media by using the well developed theory of free propagation, 
avoiding complicated mathematical calculations, (iii) The well developed parameters of free 
propagation could be directly transplanted to the SNN case to characterize the propagation 
properties, (iv) In experiments, beams are usually transmitted from free space into SNN 
media. A direct correspondence between free propagation and propagation in SNN media 
would be of practical interest for designing of experiments. 

In this letter, we describe our research into the relation between the SNN propagation and 



the free propagation, and describe the propagation in SNN media using the well-developed 
theory and clear physical pictures of free propagation. The relation between propagation in 
SNN media and that in free space is studied using the variable transformation technique. 
This shows that the governing equation, beam solutions, and propagation properties for free 
propagation can be directly transplanted with a one-to-one correspondence to the case of 
propagation in SNN media. On the basis of the one-to-one correspondence and the Huygens- 
Fresnel integral, we develop an efficient numerical method to describe SNN propagation. We 
describe the existence conditions and possible structures of solitons and breathers in SNN 
media in a unified manner by comparing the propagation properties in SNN media with 
those in free space. As an example, the theoretical predictions are illustrated for the case of 
EHG beams. 

2. One-to-one correspondence 

We begin by using the technique of variable transformation to connect propagation in SNN 
media to free propagation. This technique is frequently used to study similaritons in nonlinear 
wave guides (e.g., |19l|20]). The variable transformation technique allows the governing 
equation [e.g., the inhomogeneous nonlinear Schrodinger equation (NLSE)] to be reduced 
to a simpler mathematical form (e.g., the standard NLSE), and the solution of the former 
can be obtained from the solution of the latter by exploiting a one-to-one correspondence. 
Because the structures of the HG, LG and IG solitons in SNN media are also the modes in 
free space [18], it is reasonable to expect that a transformation exists that can reduce the 
governing equation of SNN propagation to that of free propagation. 

In the laboratory reference frame, beam propagation in nonlocal nonlinear media is gov- 
erned by the nonlocal NLSE 

2iknQd^A + rao(<9^x + dyy)^ + 2k'^/\nA = 0, (1) 

where k represents the wave number in the media with a linear refractive index uq, An = 
71-2 / -R(r — r(j)|$pd^r(i represents the nonlinear perturbation of refractive index, where ?t,2 is 
the nonlinear index coefficient, R is the normalized symmetric real spatial response function 
of the media, and r = (x, y) represents the 2-dimensional transverse coordinate vector. 
In the case of SNN media the nonlocal NLSE can be simplified to the modified SMM [10] 



2ikd^,^ + ((9^v + dy>y>)^ - k'^-f^Por'^^ = (2) 

in a new reference frame that moves with the center of mass, 

z' = z, r' = r — rc{z). (3) 



The adoption of this reference frame is important for the input fields whose transverse spatial 
momentum is nonzero [10]. In this reference frame the field becomes 



^(r', t!) = A{r' + Tc, z') exp[- 



ikM ■ (r' + r. 



2^/ 



+ 



ikM^z 



(4) 



In Eqs. (l2])-(j4j), 7 is a material constant, Pq = J \^\'^d'^r' is the input power, rc{z') = rc(0) + 
Mz'/Po is the center of mass of the beam, M = {i/2k) J {AV±A* - A*V±A)dxdj is the 
transverse spatial momentum, r' = {x',y'), and Fc = {xcVc)- 

To connect the governing equation of SNN propagation with that of free propagation, we 
adopt the transformations 



z' = 2;co[arctan( J-) + an] 



(5) 



l*(s,C) 



where w,{C) = w,o[l + {C/zco)^]'/\ RdC) = C[l + (^co/C)'], ^co = kwl,, w,o = {e^^Po)-'/\ 
a = 0, 1, — 1, 2, —2, ■ ■ -, and s = (/i, u). Then Eq. ([2]) reduces to 



{^^^. + d,,)^ + 2ikd^^ = 0. 



(6) 



Equation ([6]) is the well-known paraxial diffraction equation that governs the paraxial prop- 
agation of monochromatic beams in free space. Thus, Eqs. ([2])- ([6]) establish the one-to-one 
correspondence between the beam solution in SNN media and that in free space: 



$(r',^') = ^i^2X^(Fir',F3 



(7) 



where 



-ir[l+tan2(^)]i 
exp{ ^^(^^)^ 



22,o[tan(f-)+l/tan(:;£-)]- 



Fiiz') = 
F2{r',z') 

Fsiz') = 

«(^') = H£-arctan[tan(£)]} 

^co(-fo) 



2cotan(f- 



2cO 
1 



P07 

Note that z^o is not a constant, but varies with the input power Pq. Equation ([7]) connects 
the propagation in SNN media with free propagation. The numerous free space monochro- 
matic beam solutions and propagation properties can be conveniently transplanted to their 
counterparts in SNN media using Eq. (JTj). 

Since a one-to-one correspondence exists between the beam solutions in SNN media and 
those in free space, a general comparison between the propagation properties in SNN media 
and those in free space would be constructive. We therefore compare three beam propagation 
properties. 
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Beginning with beam patterns, we see from Eq. ([7]) that the beam in SNN media evolves 
periodically with a period Az = 2ttZco. For convenience of discussion, we divide each period 
(from z' = (2a — l/2)7iZco to z' = (2a + 3/2)7rZco) into two half-period. In the leading 
half-period (from z' = (2a — l/2)7iZco to z' = (2a + 1/2)7tZco), the evolution of the pattern 
is a condensed configuration of that in free space from — oo to +oo. There is a one-to-one 
correspondence between patterns in SNN media and those in free space, i.e., the pattern shape 
at the cross section z' in SNN media is the same as that at the cross section ( = Zco taii{z' / Zco) 
in free space. At special cross sections where z' /zco — 2a7r = — 7r/2, — 7r/4, 0, 7r/4, 7r/2, the 
beam pattern shapes in SNN media are respectively the same as that at C = — C)0, —Zco, 0, Zco, 
-l-cxD in free space. In the trailing half-period (from z' = {2a+l/2)7TZco to z' = {2a+3/2)7iZco), 
the beam patterns are the reverse of those in the leading half-period, and the evolution of 
the pattern is corresponding to that of a inverse field ( \l/(— s,^)) in free space. In fact, 
patterns of most beams are symmetrical. For these beams, the beam patterns in the trailing 
half-period are the same as those in the leading half-period, therefore the period of pattern 
evolution reduces to Az = nZco- 

The second beam property we compare is the beam width. Although the pattern shape at 
the cross section z' in SNN media is the same as that at the cross section ( = z^o taii{z' / z^o) 
in free space, the beam width in SNN media decreases by a factor of |Fi(2')| compared with 
that in free space. Explicitly, 

V J\^-z \Fi{z')\ ^ ^ 

where the superscripts (f) and (s) refer to the cases of propagation in free space and in 
SNN media, respectively, and w^^^ and w^^^ are the corresponding beam widths. Equation 
dH]) results from the self-focusing effect of SNN media. Accordingly the amplitude increases 
by a factor of |-Fi(2')|, in accordance with energy conservation. 

The final propagation property we compare involves the cophasal surfaces. For con- 
venience of discussion, we assume the radius of the cophasal surface at the cross section 
( = Zco tan(z' / Zco) in free space is R^'^^s, C), so that in free space the phase variation across 
the transverse plane can be written as exp[ikr'^ /2R^^^]. In SNN media, according to Eq. ([7]), 
the phase variation across the transverse plane at the corresponding cross section z' = z' 
would be exp[ikr'^ /2R^'^\z')], where R^^\z') is the radius of the cophasal surface in SNN 
media, and is given by 

R^'\'') = F.i,)^ F.i.r (10) 

R(/)(s,C) Rc(0 

In the special case where R^^\s^ Q = -Rc(C); R^^\z') approaches infinity, the cophasal surface 
remains planar upon propagation. 



The evolution of the beam width and of the cophasal surfaces is interdependent for SNN 
propagation. For example, when a HG beam is input at the waist and the relation zr = Zco 
(where zji is the Rayleigh distance) is satisfied, the free propagation increases the beam 
width by a factor of Fi(z') (i.e., w^^\C) = Fi{z')w^^\0)) , so that in SNN media the beam 
width remains invariant during propagation (i.e., w^^\z') = w^^^ Fi{z') = w^^\0)). Simulta- 
neously, the cophasal surface remains planar during propagation, because R^^^^s, Q = RdC) 
is satisfied. If the beam is not input at the waist and/or the relation zji = Zco is not satisfied, 
the beam width and the cophasal surface both evolve periodically during propagation. This 
property is important for the existence of solitons and breathers, as will be discussed below. 

3. Numerical description of SNN propagation 

In real optical systems, there exist many beams with irregular amplitude and phase profiles. 
Numerical methods play an important role in studying these situations. The most popular 
numerical method is the split-step Fourier method (SSFM). The SSFM is based on dividing 
the propagation length into a large number of segments and assuming that in each segment 
the diffractive and nonlinear effects are independent and may be calculated separately. Al- 
though the SSFM is much faster than the finite-difference method, it is still time-consuming, 
because of the large number of segments which is required to ensure the accuracy. Since in 
the SNN case the nonlocal NLSE can be simplified to the SMM model ( connected to free 
propagation), we develop a simple numerical method that is more efficient than the SSFM. 
In the case of free propagation, an efficient numerical method exists based on the integral 
solution of Eq. (jH]), i.e., the Huygens-Fresnel integral [18] 

^(s, = ^ / ^(so, 0) exp[^|s - sopld^so. (11) 

Because Eq. 01 ip represents a convolution of the input field with a spherical wave function, 
the field at any plane can be obtained easily from the input plane by using the fast Fourier 
transform algorithm. 

Since there is a one-to-one correspondence between the beam solution in SNN media and 
that in free space, we get the integral solution in SNN media using Eqs. (^^ and (TTT1) . as 
follows: 

^r\z') = J^{r',r',Mr',,0)dX, (12) 

where 



V9(r ,r; 



27rw^sin(^' 



ir'"^ + ir'n - 2ir' ■ r'n sec(^) 
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Equation (fT2|) is equivalent to $(r', 2;') = Fa{$(rQ, 0)}e~*", where Fa represents the fractional 
Fourier transform of order a = z' /zco- For this reason, in a separate publication J2T] we refer 
to the propagation in SNN media as the self-induced fractional Fourier transform. 

The numerical simulation of the SNN propagation in the laboratory reference frame can 
therefore be accomplished in three steps: First, we calculate the initial center of mass rc(0) as 
well as the transverse momentum M and use Eq. (jl]) to transform the field in the laboratory 
reference frame (i.e. A{r,z)) to the field in the reference frame z' = z,r' = r — r^ (i.e. 
$(r', z')). Next, we propagate the field in the reference frame z' = z,r' = r — Tc from the 
input plane 2;' = to the later plane z' = z' using Eq. (fT2l) . Finally, we transform the field at 
the plane z' = z' in the reference frame z' = z,r' = r — Tc to that in the laboratory reference 
frame using the inverse transformation 

ikM-r ikM'^z, , , 

A(r,2;) = <l>(r-re,2;)exp[^ — ^]. (14) 

This approach provides a straightforward way to numerically propagate any field from 
the input plane to an arbitrary later plane in SNN media. Since the propagation distance is 
not required to be divided into a large number of segments, only a single fractional Fourier 
transform is required, making this approach much more efficient than the SSFM for the SNN 
case. We believe the fast fractional Fourier transform algorithm which has been developed in 
recent years (e.g. [22] and references therein) would make this approach even more efficient. 

4. Breathers and solitons in SNN media 

A special feature of SNN media is that the nonlocality supports (2+l)D solitons and 
breathers [ll[3l|6l[9] , preventing the catastrophic collapse that occur in local nonlinear media. 
Here, using Eq. ([7]) and comparing propagation properties in SNN media to those in free 
space, the existence conditions of breathers and solitons in SNN media can be described con- 
veniently in a unified manner. If a beam keeps its shape during free propagation, it would 
be a breather in SNN media, because the beam shape would remain invariant and the beam 
width as well as the cophasal surface would evolve periodically with the period l\z = irzco in 
SNN media. Furthermore, if the input power and the entrance plane are designed appropri- 
ately so that the beam width and the beam shape simultaneously remain invariant during 
SNN propagation, the breather would reduce to a soliton. 

This explains why the HG, LG, and IG solitons exist in the SNN media. They retain the 
beam shape in free space, therefore they generally evolve as breathers in SNN media. In the 
special case that the field is input at the beam waist and the input power Pq equals the 
critical power Pc (so that Zco = zji, where Pc = l/(/c^7^Wo), and z^ is the Rayleigh distance 
of the input field), then diffraction increases and self- focusing decreases the beam width by 
the same factor of ^1(2:'), and the deforming of the cophasal surfaces caused by self- focusing 
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exactly balances that caused by diffraction. Therefore the beam width in addition to the 
beam shape remains invariant, and the breather reduces to a soliton. 

Furthermore, based on the analysis above we can extend the range of breathers and solitons 
in SNN media to the input fields that are linear superpositions of the degenerate solutions 
of HG ,LG, and IG beams with the same Rayleigh distance, beam waist location, and Gouy 
phase shift in free space. These superposed fields are shape-invariant in free propagation, 
thus they propagate as solitons in SNN media when the entrance plane is located at the 
beam waist and Pq = P^, otherwise they propagate as breathers. The various structures of 
these superposed fields would greatly enrich the family of solitons and breathers. 

5. Example 

To illustrate the predictive capacity of our calculations, we take the EHG beams as an 
example. In free space, the field of the {m, n) mode EHG beam can be written as [18] 



xH^{^)u)eM-c{Os\ (15) 

where c(C) = —ik/2q{Q, q{() = (—izR, zr = kw^ is the Rayleigh distance, and the coefficient 
ip is determined by the input power through Po = J |\l/pd^s. In free propagation, the {m,n) 
mode EHG beam is shape-invariant when m,n < 1, otherwise it is shape- variant. 

Obtaining the field of EHG beams in SNN media using SMM directly is mathematically 
complicated. However, using Eq. ((Tj) it is easily obtained as: 



<D^„(r',z') = Fi{z')F2{r',z')ij[ 



Qo 



'-+1 



Q{z') 



X H^,[^C{z')F^{z')x']H^[^JC{z')F^{z')y'] 

X exp{-C(^')[i^i(^')^f}, (16) 

where C{z') = -ik/2Q{z'), Q{z') = Zcotan{z' /zco) - izR, z^o = 1/7a/Po- 

Figure 1 shows the evolution of the pattern of the (2, 0) mode EHG beam in a SNN medium 
under varying input power compared with the same situation for free propagation. We find 
that if the pattern varies upon propagation in free space, it does so in the SNN medium 
as well. Because the transverse pattern is distributed symmetrically, the pattern evolves 
periodically with the period Az' = nzco- Furthermore, at the cross section z' = {a + l/2)7iZco, 
the pattern is the same as that at the far field in free space and does not vary with the input 
power, and the beam width w is inversely proportional to y/Po, because the field at the 
cross section may be regarded as the conventional Fourier transform of the input field. At 
an arbitrary cross section where z' ^ airZco/^, the pattern shape is the same as that at the 
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Fig. 1. (Colour online) Evolution of the pattern of the (2, 0) mode EHG beam in free 
space (row 1) and in SNN media (rows 2-4). The input power of the SNN propagation is 
Pq = 0.5Pc (row 2), Pc (row 3), and 2Pc (row 4), respectively, where Pc = l/(A;^7^w;o)- ^^ 
row 1 the transverse dimension is scaled by a factor of 1/[1 + ((^/z^j)^]^/^, whereas in other 
rows it is not scaled. 
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Fig. 2. (Colour online) (a) The same as Fig. 1 except that the input field is a (1,1) mode 
EHG beam, (b) Evolution of the beam width (solid lines) and the radius of cophasal surfaces 
(dashed lines). 

free space cross section ( = Zcota.n{z' /zco)- Since ( varies with the input power, the pattern 
is different for different input power. For example, at z' = 7rZco/4:, the pattern for the input 
power Pq = 0.5Pc, Pc, 2Pc is the same as the free space pattern at C = V^zr, zji, Zf(/\/2, 
respectively. 

Because the (1, 1) mode EHG beam is shape-invariant in free space, its pattern remains 
invariant in SNN media (Fig. 2). Generally, the (1,1) mode EHG in SNN media propagates 
as a breather, i.e., the beam width as well as the radius of the cophasal surfaces varies 
periodically with the period Az = vrzco- At special cross sections where z' /zco = an /2 the 
width is either maximum or minimum (in fact, the width reaches its maximum and minimum 
alternately during propagation) , and the radius of the cophasal surfaces approaches infinity 
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(i.e. tan~^(_R) = vr/2). For the special case where Pq = Pc the relation Zco = zr is ensured and 
diffraction is exactly balanced by self-focusing. Therefore the beam width and the radius of 
the cophasal surfaces are invariant during propagation, and the breather reduces to a soliton. 

6. Conclusion 

In conclusion, optical beam propagation in SNN media is connected with free propagation 
using the technique of variable transformation. The fact that the solutions as well as the 
propagation properties in free space can be mapped to their counterparts in SNN media 
through a one-to-one correspondence makes this technique useful for investigating the prop- 
agation problems in SNN media. The efficient numerical method provided in this letter is of 
interest for the investigation of beams with irregular amplitude and phase profiles in SNN 
media. The technique provides a unified description of the existence conditions and gives 
possible structures of solitons and breathers in SNN media. The various soliton and breather 
structures predicted herein would greatly enrich the family of solitons and breathers. 

Mathematically, the modified SMM is equivalent to the famous equation for the linear 
harmonic oscillator, which is widely used in many branches of physics. Therefore, the relation 
proposed in this letter can connect the free propagation not only with propagation in SNN 
media, but also with the evolution of waves in other systems in which the governing equations 
refiect the equation for the linear harmonic oscillator. The technique described in this letter 
can be readily employed in other contexts with equivalent governing equation, such as for 
quadratic graded-index media [23], or for harmonically trapped Bose-Einstein condensates 
in the noninteracting limit j24[[25]. 
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